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Abstract. We consider the asymmetric random average process which is a one- 
dimensional stochastic lattice model with nearest neighbour interaction but continuous 
and unbounded state variables. First, the explicit functional representations, so-called 
beta densities, of all local interactions leading to steady states of product measure form 
are rigorously derived. This also completes an outstanding proof given in a previous 
publication. Then, we present an alternative solution for the processes with factorized 
stationary states by using a matrix product ansatz. Due to continuous state variables 
we obtain a matrix algebra in form of a functional equation which can be solved exactly. 

PACS numbers: 05.40.-a, 02.50.-r, 45.70.Vn, 05.60.-k 
1. Introduction 

In recent years the study of stochastic systems has become an attractive and important 
research field in modern statistical physics. This is mainly based on the fact that in 
particular a lot of interdisciplinary problems are described best by probabilistic models. 
In addition, many stochastic processes represent simple nonequilibrium systems and may 
serve as a kind of toy models for the evaluation of a still outstanding nonequilibrium 
theory. 

In this work we focus on the asymmetric random average process (ARAP) 
H]. The model is defined on a lattice and equipped with a probabilistic nearest 
neighbour interaction. This is in common with most stochastic models, especially 
with the asymmetric simple exclusion process (ASEP), e.g., and references therein, 
representing somewhat like a standard model of nonequilibrium physics. However, the 
state variables of the ARAP located at the lattice sites are continuous and unbounded 
while most of the models known deal with discrete and even finite local state spaces. 

Nevertheless, the ARAP is not an artificial construction. Many physical problems 
are rather located in continuous than in discrete space. E.g., it is closely related with 
the q-model of granular media [Q. The traffic model of Krauss et al can also be 
mapped onto the ARAP Furthermore, the ARAP may show some new phenomena 
undiscovered in discrete systems so far, e.g., in a new kind of twofoldly broken 
ergodicity has been studied. 
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The results presented in this paper are structured as follows. 

In section ^ we give a definition of the ARAP in context of the so-called quantum 
formalism of stochastic processes This is useful for the application of the matrix 

product ansatz (MPA) presented in section ^. 

In sections ^ and |^ we focus on ARAPs with interactions leading to factorized 
steady states. In |TD] these interactions have been identified, however, the calculation 
was not rigorously and involved a unproven conjecture. According to this, we complete 
the outstanding proof in section ^. Furthermore, we derive an explicit functional 
representation of the interactions, so-called beta densities, and discuss the infiuence 
of finite system sizes. 

Finally, in section ^ we apply the matrix product ansatz to ARAPs with interactions 
identified in the previous section. This approach has been successfully developed 
for quantum spin chains and stochastic systems in the last decade (see section H for 
references), however, it has always been applied to systems with discrete state variables. 
Here we enhance its validity to systems with continuous state spaces. The corresponding 
algebras become compact functional equations with closed solutions. 

In the conclusion we discuss the MPA for ARAPs with interactions that do not 
lead to product measure steady states. 



2. Definition of the process 

The ARAP is defined on a ID periodic lattice with L sites. Each site i carries a non- 
negative continuous mass variable mjGlR,}. The dynamics is given in discrete time and 
every time step t — > t + 1 for each site a random number G [0, 1] is generated from 
a time-independent probability density function (j) = (p{ri), called fraction density. The 
fraction determines the amount of mass Aj = rjmj transported from site i to site i + 1. 
The transport is completely asymmetric, i.e., no mass moves in the opposite direction 
i + 1 — s> z, and we obtain 

rrii ^ {1 - ri)mi + ri^irrii^i . (1) 

These update rules correspond to a parallel dynamics and are illustrated in figure ^ 

So there are only three parameters in the ARAP: first, the system size L and 
the mass density P = x whereby M = Yli '^i represents the total mass. Due to 
mass conserving dynamics the density is fixed and accordingly p can be considered as 
a thermodynamic variable. However, the most powerful "parameter" is given by the 
fraction density (f) which allows the ARAP to be customized to a lot of interdisciplinary 
problems. In |0 some stochastic models have been rewritten in context of the ARAP 
by the use of suitable (/)-functions, e.g., the q-model from granular media or the Krauss- 
model from traffic fiow theory. 

Now the ARAP will be presented in terms of the so-called quantum formalism for 
stochastic processes, e.g., used in We consider the orthonormal state space which is 
spanned by the continuous ket basis {|ni)} with configuration vectors m = (mi, . . . , mi) 
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Figure 1. Mass representation of the ARAP with paraUel dynamics. The height of a 
mass stick corresponds to rrii. The fragments rirrii are shaded. 

and equipped with the inner product (mi |m2) = 6 (mi — m2) (here 6 represents the 
^-function). Correspondingly states of the ARAP at time t are given by 

POO 

|P(t)) = / d^mP(m,t) |m) (2) 
Jo 

whereby d^mP(m, t) is a non-negative probabihty measure or reworded: P(m, t) gives 
the probabihty density of finding the system in the configuration m at time t. Here the 
abbreviated form jj d^x = jj . . . jj d^x with / C IR has been introduced. Furthermore 
we assume P to be normahzed, i.e. d^mP{in,t) = 1. 

Our main aim is to calculate the function P which corresponds to the solution 
of the problem. In case of the parallel dynamics given above we obtain the following 
t t + 1 map of the basis: 

|m)^T|m)= [ d-^r(/)(r) |T(r)m) (3) 
Jo 

with (f){r) = Yl - 4>{ri) and 

Tij(r) = (1 - ri)5i^j + ri„i(5i_ij . (4) 

Here T(r) represents an L x L matrix with diagonal elements 1 — rj and lower band 
entries r^. Based on periodic boundary conditions the top right entry is also unequal 
to zero, i.e. Ti_i(r) = r/,. So m ^ r(r)m is nothing else than the compact matrix 
formulation of (0). Note that T(r) operates on the space of configuration vectors m, 
whereas T operates on the state space spanned by |m). 

Evaluating P(m',t + 1) = (m'|P(t + l)) = (m'|TP(t)) by using the relations 
(0),(^ and (^) we finally obtain the master equation 

1*00 I* 1 

P(m',t + 1)= / dV / dV 0(r)5(m' - r(r)m)P(m,t) (5) 
Jo Jo 

representing the fundamental dynamical equation of the ARAP. Here the expression 
0(r)5(m' — T(r)m) represents the transition probability density from state |m) into 
state |m'). 



r2m2 



^ 13 m3 



t — t+1 
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From now on we would like to focus on steady state dynamics only, i.e., we look 
for time-independent solutions P(m) of (|^). This simplifies the master equation to the 
eigenvalue problem |P) = T|P). For ergodic ARAPs, defined by dynamics that allow 
the system to evolve from any given initial state to any final state in a finite time, the 
steady state is equal to the infinite time limit, so |P) = limt^oo |P(^))- 

Although we do not focus on ARAPs with discrete masses or other kind of 
updates in this paper, we would like to mention that it is straightforward to deduce 
the corresponding master equations from (|^). E.g., the discrete ARAP is embedded 
canonically by using ^-functions, whereas ARAPs with continuous time dynamics can 
be treated by a suitable chosen fraction density that interpolates between parallel and 
random sequential (=continuous time) updates. Both cases are explicitly treated in f^. 

Furthermore, we do not take into account systems with state-dependent fraction 
densities, e.g., given by 0(r, m) (here the probability of mass shifts also depends on the 
height of the column). Although these ARAPs may show interesting phenomena |^, an 
analytical treatment is more difficult in general. 

So, the ARAP with continuous state variables, discrete parallel update and 
state-independent fraction density (j){r) spans an appropriate framework for analytical 
research. We hope that a lot of results can be transferred to related ARAPs in a next 
step. 



3. Product measure solutions 



In this section the complete set of ARAPs with product measure solutions 
P(m) = n^( 



irrii 



(6) 



is presented, i.e., we determine rigorously the set A4 of all 0-functions leading to 
factorized mass distributions. Basically, we complete the proof of [|l^ which includes a 
conjecture based on high order computations so far. Furthermore, the fraction densities 
of Ai are given in a closed form in contrast to our previous paper [10| where the (l>- 
functions are in moment representation only. Finally, we briefly discuss the influence of 
the system size L. 



3.1. Explicit form of M. 

In [|10[ it is shown by a constructive approach that is a two parametric set of 0- 
functions, determined by the first and second 0-moments /Ui and /i2 which are defined 
by //n = /q*^ dr (/)(r)r". These free parameters and /i2 have to be chosen with respect 
to the general moment properties 

1 > ;Ul > ^2 ^ (7) 

only. Then the higher moments are uniquely determined by and ^2 through 

_ r(n + A) r(A) 
r(A) r(n + A) 
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with 

A = fii— — ^ and A = — — ^ . (9) 

From now on we neglect for simplicity the special case /X2 = /^i corresponding to 
0(r) = (5(r— /ii) and leading to product measure solutions generated by P{m) = 5{m—p). 
In particular, these singular 5-densities will not be considered as elements of M.. For 
completeness we present the factorized mass densities associated with (||), also deduced 



m 



A^ 1 fmV'^ _y_ 



Please note that the uniqueness of the higher moments fin>3 has already been proven 
p!0[] . Correspondingly, the densities given above span the complete set Ai. However, 



m 

the explicit form (|[) is a conjecture that has only been shown rigorously for n < 10. 
According to this, our approach in this section is canonical: We first derive explicitly 
the fraction densities related to (H) and after that we show that these 0-functions really 
lead to factorized mass distributions. This will complete the outstanding proof. 
We start by rewriting (§) as a recurrence relation: 

fJ-n+l = ^ ^ fJ'n with yUo = 1 • (H) 

n + X 

Then we define the generating function 

^(^) - E ^/^«^" • (12) 

n 

From pn+i ^ /Un we derive that F is an entire function and consequently well defined. 
In addition, the relation 

F{s) = [ dr0(r)e''' (13) 
Jo 

holds and F{is) is the characteristic function of (p. 
From (0) and (|12D the differential equation 

sF"{s) + (A - s)F'{s) - ~XF{s) = (14) 



is derived after some algebra which is nothing else than Kummer's equation 11|. This 



ordinary differential equation of second order is elaborated very well, e.g., it appears in 



context of the hydrogen atom, and so we may rely on a huge pool of known results |11 



For linear differential equations of second order there are always two independent 
solutions. But here, only one of them, the so-called Kummer M-function M(A, A,s), is 
analytical in s = as long as A is not a negative integer. This is ensured by (0) and (^. 

Consulting [|1^] yields for the special case Re A > Re A, which is also satisfied here, 
the integral representation 

M(A, A, s) = f dr ^^}^^ r^-^(l - rf-'-'e^^ (15) 
Jo 1 (Aji (A — Aj 
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which finally gives an explicit representation of by comparison with (|T3|). 

At last, this functional solution is rewritten in terms of parameters a and b instead 
of A and A to simplify the representation of A4. By the transformation a = A and 
a + b = X a symmetric form of the fraction densities is achieved and we get 

^ = (0a,6(r) = — ^r"'^(l-r)^-ia,6GR+| . (16) 
I B{a,b} J 

These are so-called beta densities which are very common in probability theory |T^. At 
this, the normalization constant is given by the beta function 



Although the mean field models (piq ) are parameterized by a two-dimensional 
manifold, the associated mass distributions (|10D are connected to a one-dimensional 
parameter space only, i.e., several <pa,b functions yield identical stationary states with 
X = a + b. However, we have to keep in mind that this conclusion concerns only the 
steady state. The relaxation into the stationary state could differ completely. 



Now we briefly discuss the product measure ARAPs (|T^) into more detail. In 
dependence of the fraction density the process may behave from even critical (a + 6 — > 0) 
to deterministic (a -|- 6 — > oo). Here critical ARAPs are characterized by an algebraic 
mass decay [Q. They are realized by dynamics that either shift the total mass located 
on a lattice site or forbid the transport. 

In general, there are three classes of A^-densities: continuous functions with 
0(0), 0(1) G {0,1}, single peak functions with 0(0) = oo or 0(1) = oo and double 
peak functions with 0(0) = 0(1) = oo. Only double peak functions may result into 
mass densities P(m) that diverge for m — ^ 0. This reflects an almost critical behaviour 
because either transports of no mass (r = 0) or the total mass (r = 1) are preferred. 
However, all classes may generate Gaussian-like mass densities with an algebraic increase 
for small masses. Furthermore, it is interesting that continuous and single peak densities, 
resp. single and double peak densities, may lead to identical mass distributions, whereas, 
this is impossible in case of continuous and double peak ARAPs. 

For a = 6 = 1 we obtain the simplest version of the ARAP with uniform distribution 
0(r) = 1. We will refer to this system as the free ARAP. In particular, there is no explicit 
truncation which can forbid or suppress a transport of fractions bigger than a critical 
value, e.g. discussed in 0. 

3.2. Completion of the proof 



Now we prove that the densities (16) really lead to factorized mass distributions 



according to (|ToD- In 1131 ^ simple criterion has been presented which makes it possible 
to determine and verify mean fleld solutions: Only product measure mass densities 
satisfy the equation 

Fq(si,S2)=Fq(si,0)-Fq(0,S2) (18) 
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with 

Fq(si, S2) = [dr (P{r) Q ((1 - r) Si + rss) (19) 
Jo 

and single site Laplace-transform Q{s) = dmP(m)e~™"^. Hence, (|18D represents a 
sufficient criterion for testing the validity of product measure solutions. In the following 
paragraphs we will show that the mass densities (|T0|) meet the condition ([T8|). 



The Laplace-transform of (|T0|) is given by 



p x-A 



Q{s) =(1 + ^s) . (20) 

Then, using (|16D and (1^) with (ID leads to 

, T(a + b) ^a-i/^ _ 
«^ io r(a)r(6)(i + £[(i_^),^ + ^,2])^ 

Evaluation of the r.h.s. of (pi]) can be done by use of the so-called Feynman parameters 
which are a well-known tool in field theory. They are nothing else than the formula 

En 

Here a, are real and positive, whereas Xi may be complex. A derivation of (0) is given 
in 



13 1 . Now we simply apply the values 



n = 2 , Xi = 1 + , ai = b and 02 = a , (23) 
A 

integrate over ri, relabel r2 — r, and finally equation (P^) rereads as 



p \-b / p \-« r\ T(a + b) ra-ui_ \b-i 
1 + tSi) (I + TS2) = / dr j Zn\ ^ r- 24 

AV V aV J, r(a)r(6)(i + £[(i_^),^ + ^,j)^ ^ ' 



Together with (|2l|) we derive 



Fq(si,S2)= (l + ^Si)"'(l + ^S2) . (25) 



From this follows directly the validity of condition (|TH|) which completes the proof. 

The formula ( p2D has already successfully been used in context of the q-model |^[T^ 
As mentioned in the introduction this fundamental process of granular media is strongly 
related with the ARAP. However, only symmetric ARAPs defined by symmetric fraction 
densities 0(r) = 0(1 — r) can be mapped onto the q-model |1|,|§1. Correspondingly, the 
application of Feynman parameters has been generalized to antisymmetric 0-functions 
here. 

3.3. Finite systems 



The calculations in |TD| have only been focused on systems in the thermodynamic limit 
L ^ 00. But, also in case of finite systems the mean-field criterion (|T^ is valid and 



sufficient which can be shown in almost the same manner as done in 10 . So ARAPs 



with beta densities are of product measure form irrespective of the system size L. 
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However, for L < oo the configuration space is restricted to the hyperplane 



L 



Fl(M) = |m I ^ mi = m| (26) 

due to mass conserving dynamics. Accordingly, the mass density has to be renormalized 
and we obtain 

^(m) = I ^ nf=i ^A(mi) for m G Fl{pL) ^^7) 
1 dsG 



with 



Z = Z(A,p,L,pL)= I d^mnPAK) = -^^^e-^^ (28) 

Thus, projection onto the F^^pL) surface, i.e. fixing the total mass, corresponds 
to shifting our focus from grand-canonical to a canonical point of view whereby Z 
corresponds to the canonical partition sum. A detailed calculation of Z is given in 
[Appendix A| . 

One should keep in mind that the exact solutions (|27D are still of product measure 
form if restricted to Fl[pL). However, this coincidence with the L = 00 case is only 
formal, e.g., the one-site mass density is not simply given by Z^t: P{m). One has to 
take into account the additional interaction induced by the restriction m G Fi^pL). So 
we derive by using the relation ( [A.l| ) again: 



K) = / d-^m ip.(m.) n F.(m,) = ^^^'^'(;^^;'^'^'ip^^ ^A(m.) .(29) 

FL-i{pL-mi) 

4. Matrix product ansatz 

In this section we solve ARAPs with 7W-densities by using a matrix product ansatz 
(MPA) . This technique has been initially introduced for calculating exact ground states 
of quantum spin chains [|T5|. Shortly after, Derrida and coworkers have successfully 
applied the MPA to a nonequilibrium system, namely the ASEP with random sequential 
dynamics |]TB|]. Meanwhile the MPA has been evolved to a standard tool for one- 
dimensional stochastic models, e.g. ^ |T^, and references therein. However, 



its field of application so far is mainly restricted to variants of the ASEP (different 
updates, local defects or two species of particles). 

In general the MPA is applied to systems defined on a two-dimensional local state 
space, i.e., sites can be vacant or occupied. In this case, manageable sets of algebraic 
objects (corresponding to the local states) and algebraic relations (refiecting the local 
dynamics) are obtained. For example the ASEP provides one condition, pDE = D + E, 
where the objects E and D correspond to holes and particles. Nevertheless, it is rather 
complicated to find (matrix) representations of that poorly defined algebra 0]. An 



Matrix product approach for the asymmetric random average process 



9 



extension to a model with an arbitrary, but still finite, number of local states is given 
in 



20 



Here we apply the MPA to a stochastic system with continuous state variables. As 
seen in the previous section, mean field is exact for fraction densities taken from Ai, 
i.e., the corresponding algebras have one-dimensional representations and are given by 
a functional equation. 

We derive mass solutions in agreement with the results given in section So an 
alternative approach for the calculation of steady states is presented in this paper. This 
may show new perspectives in the treatment of other ARAPs which are unsolved so 
far. E.g., the free ARAP with continuous time update still lacks an exact description 
of the steady state. Also the state-dependent models given by mass-dependent fraction 
densities 0(r, m) could be treated by the MPA. 

Finally, our work extends the scope of the MPA to systems with unbounded and 
continuous state spaces! Thus, we have to deal with functions and functional equations 
instead of discrete objects. For this, we may fall back on the well elaborated field of 
functional theory. So, although things get in principal more complex - measured in 
degrees of freedom - life becomes easier. 

4.I. Continuous algebra 

In this section we derive the matrix algebra of the ARAP in the thermodynamic limit. 

We start with the "defect" matrix product ansatz for backward sequential dynamics 
T7| . A definition of the backward sequential update is presented in the [Appendix B| 



whereas a detailed introduction to the MPA for stochastic systems can be found in 
the references given in the introduction of this section. Therefore, we assume the local 
interaction to obey 

t{A^A)=A®A (30) 

with 

/•oo /"OO 

A = dmD(m) \m) and A = / dmD(m) \m) . (31) 
Jo Jo 

Here \m) spans the infinite and continuous local state space of a single site and the 

tensor product is defined as usual, i.e. 

/•oo roo 

A0A= / dm / dm D{m)D{m)\m,m) , (32) 
Jo Jo 

whereby we have used |m,m) = \m) ® This is in accordance with the notation 
introduced in section ||. Note that the algebraic objects D and D depend on a continuous 
parameter m reflecting the mass located on a lattice site. 
It is easy to see that 



P)^ = tT\ A0A0...®A0A\ (33) 
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represents a steady state under backward sequential dynamics 

%s = ^L, 1^1, 2^2,3 • • • ^L-l,L • (34) 

The operator is defined by the local interaction t acting on sites i and i + 1. For 



further information about update procedures please refer to [0 . 

Note that trace operator tr and time evolution operator Tbs commute because tr 
acts on the auxiliary space of the algebraic objects whereas Tbs is defined on the state 
space. Here the trace operator shall ensure translational invariance of the steady state. 

In the thermodynamic limit the parallel update corresponds to the backward 
sequential update (see [Appendix B| ), so the defect can be neglected and we obtain 

|P) = lim |P)^ = tr (A® A® . . .) . (35) 

L— >oo 

Now we give the explicit definition of the local dynamics. First, we change notation 
(or, more formal, the basis) uniquely in the following way 

\m,m) — » \m + m,m) (36) 

because the local dynamics is mass conserving. Then, the representation of the local 
interaction reads as follows: 

t\s,m)= / dA/(A,m) |s,m - A) . (37) 
Jo 

Here / represents the so-called fragment density, giving the probability (density) that a 
fragment of size A is broken of a stick with mass m. The densities / and <!) are simply 
related by 

f{A,m) = -c^(-) . (38) 

Using (|30D -(^) with (|36|)-(p8D yields after some calculation the matrix algebra in 
form of a functional equation: 

/ dA -01 -]D{s- A)D{A) = D{m)Dis-m) . (39) 

Jo s-A \ s-AJ 



For a detailed derivation we refer to [Appendix C| . Please note that we are confronted 



with one condition and two objects only. This is remarkable due to the immense degrees 
of freedom. In discrete systems every local state corresponds to a single algebraic 
object, e.g., in the ASEP we deal with an object E for empty sites and D for occupied 
sites. Accordingly, unbounded local state spaces yield an infinite number of algebraic 
objects. Additionally, the number of conditions derived from the local dynamics and 
defining relations between the algebraic objects depends quadratically on the number of 
local states. So, increasing the degrees of local freedom makes it more difficult to find 
representations of the algebraic objects that fulfill all conditions. Even more fascinating 
is that in case of continuous system things become easier again! The algebraic objects 
condense to functions and many conditions summarize in one functional equation. 
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4-2. Constructive solution for the special case 



In general D{m) and D{m) are arbitrary algebraic objects and sometimes it is 
possible to derive information about the underlying system without finding a concrete 



representation, e.g., this has been done for the ASEP in ||T^ where some quantities 
like the flux have been calculated by recurrently solving the algebraic relations. 
Usually one tries to find a matrix representation fulfilling the dynamical conditions. 
Here one distinguishes one-dimensional and higher dimensional representations because 
one-dimensional representations correspond to product measure solutions (||) always, 
whereas higher dimensional solutions incorporate additional correlations and show that 
the process cannot be solved by a mean field ansatz. 

In this subsection we will focus on ARAPs with ^i_b-functions, i.e. a = 1, given by 

0i,,(r) = 6(1 - rf-' (40) 

and part of the class Ai. Without any further input from section |^ we will derive the 
according mass densities by the matrix product technique, i.e., by solving the algebraic 
equation (|39|) . Since we are looking for product measure solutions, we assume a one- 
dimensional representation of D and D. Even without knowing that such solutions 
exist, one would always try to solve a problem by the simplest ansatz. If this approach 
fails, we would look for higher dimensional representations. 

This means that D and D are nothing else than functions with one real parameter. 
In particular, D and D commute and we assume them to be differentiable. 

Using 01 with (|39D we obtain 



ps—m 

hi dA{s-A)-^D{s-A)D{A) = m^-^D{m)D{s-m). (41) 
Jo 

Then differentiating with respect to m generates after rearrangement 
(1 - b)D{m)D{s - m) + bD{s - m)D{m) 

= mD{m)D'{s — m) — mD'{m)D{s — m) , (42) 

a differential equation with two functions and two variables s and m. For s = 2m the 
explicit 6-dependence vanishes and (H2|) reduces to 



_d_ 

dm 



■ D(m) 
In — r 
Dm) 



1 

m 



(43) 

uyiiLj \ in 
Directly the relation 

D{m) = CmD{m) (44) 

is obtained. Inserting this solution in (^) enables us to extract D from the algebra and 
after some calculus we derive 

h{m) - h{s -m) = {b-l) (45) 

m[s — m) 



with 



h{z) = -^\nD{z) . (46) 
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For the free ARAP, i.e. 6 = 1, the r.h.s. of p5| ) is equal to zero. Then h has to 
be constant because has to be valid for all ^ m ^ s < oo. Thus, we get 
D{m) = Cexp{iim). For arbitrary b we rewrite the r.h.s. of (|45| ) by expansion into 
partial fractions and achieve the difference equation 

h(m)-h(s-m) = (b-l)(— —] . (47) 

\m s — m J 

Its general and unique solution is given by h{z) = fi + {b — l)z^^ (at this we used 
the same argument as for the 6=1 case). By the help of definition (^) we obtain 
D{m) = Cm}'~^ exp(/im) and finally D{ni) = Cmf' exp{^m) by (0). 

In case of a one- dimensional solution the trace operator in (^) becomes redundant 
and the identity P(m) = D{rn) is valid. Using the boundary conditions dmP{m) = 1 
and Jq°° dmP{m)m = p, we determine the constants C and /i and obtain the single-site 
mass density 



:i + bY+^ 1 f 



m 



P{m) = ^py^- ( - ) e-(^+^)^ . (4^ 

r(i + 6) p VP 



which is in perfect accordance with ([ToD . 

We emphasize that the matrix product ansatz offers a big advantage compared 
to other approaches dealing with product measure solutions: we do not have to prove 
the exactness of the solution (|48|). Usually one assumes the master equation to be 
solved by a solution of type (|]), i.e., we make a so-called mean-field ansatz and look 
for solutions 0, ||,[10|. The next step is to show that the mean-field ansatz is really 
correct, i.e., we have to prove that all higher correlations or joint probability densities 
decompose. Although in |T0[ this extensive task has been reduced by deduction of the 
criterion (|1^), the step of testing remains: mean-field solutions have to satisfy ( [I8| ) for 
all si and S2. 

However, in case of the matrix product technique the proof of exactness is delivered 
for free. The construction of the steady state ensures its exactness (see subsection ^T| ) 
and the dimensionality of the representation gives information about the correlations. 



4-3. Approach for the general case (j)a,b 

Here the matrix algebra for beta densities with arbitrary b and a G N is considered. 
As in the previous section we derive the corresponding differential equations involving 
D and D, however, these become more complex and we have to treat them for each a 
separately. We present the a = 2 case in detail and refer to the problems in finding 
closed solutions for arbitrary a. 

Starting point is (|39|) with fraction densities ([T6|) . Introducing H(z) = z^^^"'^''^ D(z) 
we rewrite the functional equation to 

ps—m 

/ dA{s-m- Ay~^H{s- A)D{A) = B{a,b)m^-^D{m)D{s-m) . (49) 
Jo 
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Now we assume a to be integer and derive the relation 

/ dA(s - m - A)''-^H{s - A)D(A) = (-)"(a - l)\H{m)D{s - m) . (50) 

For further information please consult [Appendix D| . So, by an a times differentiation 
the integral expression disappears and combination of ( ^9]) and (0) yields the equation 



ga 



- [m^-^D{m)D{s - m)] = (-)" ^^"J" ^^ m^~("+^^D(m)D(g - m) .(51) 



dm'' L V / V /J V / Y{b) 

For a = 1 we obtain (^21) . Defining G{z) = z^^^D{z), evaluating the expression on the 



m-^D{m)G{m) . (52) 



Y^(^\-fD^^\m)G^'^-^\m) 

u—n \ / 



l.h.s. of (51) and setting s = 2m results in 

_ r(a + 6) 

r(^) 

fc=0 ^ ^ ^ ^ 

Here the upper index denotes the k^^ derivative. Thus, equation ( |52D represents a 
homogeneous linear differential equation of order a in D and G respectively. In terms 
of F = {D,G) it is even nonlinear. 

Henceforth, one may apply the following strategy as done in the previous section: 
first we try to solve ( |^ in D, i.e., we obtain a relation D{m) = J^{D){m). Then, we 
insert this result into (^T]), obtain an equation with D only and try to solve it. However, 
this approach becomes more difficult for a 7^ 1 because the corresponding differential 
equations are more complex. 

As an example, we would like to discuss the case a = 2 into detail. For this (^) 
transforms into 

D(m)G"(m) - 2D'(m)G'(m) + D"(m)G(m) = D(m)G(m) . (53) 

m^ 

To simplify ( |53D we make use of the substitutions d = (InD)' and g = (InG)' so that a 
nonlinear differential equation of first order is achieved finally: 

[d{m) + g{m)] + [d{m) - g[m)f = Ml±^ . (54) 
dm m"^ 



Unfortunately, we are not able to give a general solution of (|5^ ) representing a Riccati 



equation in a = d — g [21|. This is based on the fact that g is also an unknown 



function. Therefore, we cannot apply elaborated solution schemes for this type of 



equation here |21 



Nevertheless, we are able to present a special solution of (0) given by 

d{m) = h /i and g(m) = h /i (55) 

m m 

with (k — R)"^ — {k + R) = b{l + b). Going back, exact expressions for D and D are 

derived, especially Dim) = Cm'^ exp{fim). Although these functions do not fulfill the 

overall condition ( |5TD for a = 2, this heuristic approach yields the correct form of D 

(see below), but an erroneous defect term D. 

However, we have not exploited the relations between D and D to the full. 

According to this, we are confronted with three parameters C,k and fi now, whereas in 
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the previous section only two unknown constants occurred. To avoid this situation, the 
second mass moment 



Jo ^1-^2 

is used as a third condition. It has been calculated exactly in [Q. With a final 
reassignment D P the solution 

P,„,) . (il^i f !^V"eH«)? (57) 



r(2 + 6) p\p 

is obtained which is in perfect accordance with (p^O]). 

So the case a = 2 differs completely from the case a = 1 presented in the previous 
section: it seems to be impossible to derive both, D and D from the algebraic equation 
(p9|). We have to keep in mind that the defect term D is merely an auxiliary construct 
and not part of the final solution (^). Accordingly, there could be some freedom of 
choice for the defect. E.g., a lot of matrix product solutions for stochastic systems are 
characterized by an a priori choice of the defect terms [P,|l^,[l9|. This simplifies the 



algebra and leads to the proper solutions. In this way, (|^) represents a suitable choice 
of the defect function g. 

In principle, one could try to solve the algebra for a = 3, 4, . . . explicitly and propose 
the general form (|10D for arbitrary a. But this approach is discontenting. On the one 
hand, we do not know about the complexity and resolvability for increasing integer a. 
On the other hand, we do not provide a method of resolution for not-integer a. Finally, 
we are still lacking of an enclosing MPA calculus that is valid for all Ai densities. 
However, the MPA might nevertheless be useful for special ARAPs with non factorizing 
mass densities. 

5. Summary and Outlook 

In this paper we have focused on the asymmetric random average process under parallel 
dynamics. In the first part (section ^) it has been rigorously proven that only ARAPs 
equipped with beta densities and ^-densities lead to product measure solutions. Here, 
our calculations have completed a derivation given previously in |TD]. In a second 



part (section ^) we have presented an alternative approach, a matrix product ansatz, 
for deriving the mass densities under use of beta densities. Restricted to fraction 
densities (j)i^b we have given a complete calculation of the associated mass densities. 
For the remaining densities (pa^b with a 7^ 1 we have outlined an approach based on 
differential equations. Altogether, we have extended the scope of the MPA to systems 
with continuous and unbounded state variables. 

However, due to the restriction to ARAPs in class Ai we deal with a one- 
dimensional representation of the algebraic objects D and D only. In a next step, 
we should concentrate on ARAPs with non-vanishing correlations which is associated 
with higher dimensional representations of the algebraic objects. Especially, we believe 
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that the MPA could be the appropriate tool to solve the free ARAP with (p{r) = 1 under 
continuous time dynamics. Although this process is given by a rather simple fraction 
density it does not belong to the class of mean-field models 0,0. Furthermore, ARAPs 
with state-dependent fraction densities (j){r, m) could be treated with the MPA. This 
could bring forward the analysis of truncated processes. On the other hand, we could 
ask for 0-functions that lead to solvable algebras, i.e., solvable functional equations. 

In addition, the question arises if the algebraic objects of correlated ARAPs are 
always representable by matrices, i.e. D{m) = {Dij{m)). Or is it even possible to 
construct ARAPs that are solved by algebraic objects with continuous representations? 
Finally, one could identify the class of ARAPs whose stationary states may be written 
as matrix products. This would be very similar to [jl8[ where it has been shown that 
stochastic processes with a finite number of local states, a finite range of interaction and 
continuous time dynamics are (formally) always solved by MPA. 
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Appendix A. Partition sum Z 

We present a proof of the relation 

Z(A,L,p,M)^ j d^mn^^^K) = ^r(AL) (^A^^ e'^^ . (A.l) 

Fl{M) 

The superscript in reminds of the p dependence. The following calculation uses the 
Fourier representation of the 5-function and the relation between ([T0|) and (pOD, i.e., the 
principles of Laplace transformation: 

Zi\,L,p,M)= / d^mr[F((m,)5(M- Vm,) 

1 r°° 

= 9^ dpe'^^'l[Q{{ip) 

J -oo j 

1 r°° 

So the r.h.s. of (|A.1|) is nothing else than a single site mass density (0) for the total 
mass M that is rescaled according to p ^ pL and A XL. 
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Appendix B. Backward sequential update 

In the backward sequential update one starts updating an arbitrary pair of sites (i, i + 
and applies the local update rules, 

rrii ^ {1 - ri)mi and rui+i ^ nii+i + rirrii , (B.l) 

from right to left, i.e., opposite to the direction of the mass transport, under 
consideration of periodicity until reaching the initial pair. Thus, for finite L the 
backward sequential update differs slightly from parallel dynamics because the last pair 
of sites {i + l,i + 2) to be updated holds the site i + 1 that has been updated already. 
However, in the infinite system L — * oo this "error" vanishes because the difference 
between both update procedures is local and not extensive with L. Correspondingly the 
parallel update is equivalent to the backward sequential dynamics in the thermodynamic 
limit. 



Appendix C. Derivation of the algebra 

Putting (pO|)-p^) and (p6D-(^) together we obtain the expression 

t(A®A)=f ds f dm f dA —) D(m)D(s -m)\s,m- A) . (C.l) 
Jo Jo Jo ^ 

We introduce the linear substitution x = s — m and y = m — A that comes along with 

the trivial Jacobian determinant —1. Therefore, the transformation of ( |C1| ) becomes 

rather simple: 

t(A(^A)= f ds f dy f \x ^—(p(^—^^—^] D(s-x)D(x)\s,y) . (C.2) 
Jo Jo Jo s-x \ s-x J 

Here the change of variables generates a structure according to ds Jq dy . . . \s,y). 
Since \s,y) are linearly independent we are able to reduce the condition ([30| ) to the 
functional equation (|39|). 



Appendix D. Derivation of relation (|50|) 

The essential component of the calculation is the formula 

d r r d 

-7- j dyf{x,y) = f{x,x) + dy—f{x,y) (D.l) 
dx Jo Jo dx 

which can be verified easily. 

By the help of it is straightforward to prove ( ^Ot ) by induction. In particular, 

the correspondent term of f{x,x) reduces to zero for a > 1. However, for a ^ N 
induction finally leads to the case < a < 1 which cannot be treated by ( p.l|) because 



the counterpart of f{x,x), i.e., basically (s — m — A)*^ ^ at A = s — m, is divergent here. 
Therefore, we have restricted to a of integer form. 
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